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Abstract 

The 2-point correlation form factor, A^t), for small values of r is 
computed analytically for typical examples of pseudo-integrable sys- 
tems. This is done by explicit calculation of periodic orbit contribu- 
tions in the diagonal approximation. The following cases are consid- 
ered: (i) plane billiards in the form of right triangles with one angle 
ir/n and (ii) rectangular billiards with the Aharonov-Bohm flux line. 
In the first model, using the properties of the Veech structure, it is 
shown that A2(0) = (n + e(n))/(3(n — 2)) where e(n) = for odd 
n, e(n) = 2 for even n not divisible by 3, and e(n) = 6 for even n 
divisible by 3. For completeness we also recall informally the main 
features of the Veech construction. In the second model the answer 
depends on arithmetical properties of ratios of flux line coordinates to 
the corresponding sides of the rectangle. When these ratios are non- 
commensurable irrational numbers, A2(0) = 1 — 3a + 4a 2 where a is 
the fractional part of the flux through the rectangle when < a < 1/2 
and it is symmetric with respect to the line a = 1/2 when 1/2 < a < 1. 
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The comparison of these results with numerical calculations of the 
form factor is discussed in detail. The above values of 1^(0) differ 
from all known examples of spectral statistics, thus confirming ana- 
lytically the peculiarities of statistical properties of the energy levels 
in pseudo-integrable systems. 
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1 Introduction 



The statistical properties of quantum systems attracted wide attention in the 
last years (see e.g. Jl|). The investigation of many different models had led 
to a few accepted conjectures which relate statistical distribution of quantum 
energy levels with general properties of corresponding classical motion. For 
generic systems these conjectures are the following: for chaotic systems the 
level spacing distribution follows the Random Matrix statistics [[J, 0; for 
integrable systems it follows the Poisson statistics ||. Both conjectures are 
supported by a lot of numerical evidences and by some analytical arguments 

1-0- 

These well-established conjectures are applicable only to completely chao- 
tic or integrable models. But there are systems which are neither chaotic nor 
integrable. Noticeable examples of such systems are plane polygonal billiards 
with all angles, ati, commensurable with ix 

OLi = 71 — , (1) 

Hi 

where rij are co-prime integers. In such systems all trajectories belong 
to a 2- dimensional surface of genus 

iV \ mi — 1 

i 

where N is the least common factor of the rij (see e.g. [pj|). The case where 
all Ttii — 1 corresponds to g — 1 (i.e. to a torus) which is integrable. If some 
rrii > 1 trajectories belong to a higher genus surface and, consequently, the 
system is not integrable (at least in the usual sense) but it is not chaotic either 
since all trajectories belong to a 2-dimensional surface and cannot cover a 
3-dimensional energy surface ergodically as is required for chaotic systems. 
For such reasons these systems are called pseudo-integrable. 

A natural question appears: what is the spectral statistics of pseudo- 
integrable systems? Numerical calculations |PJ| -||13| clearly demonstrated 
that statistical properties of such systems differ from standard examples but 
have many points in common with the statistics of the 3-dimensional An- 
derson model at the metal-insulator transition point [14]. The full analytical 



approach to this question meets with difficulties related mostly with the exis- 
tence of quickly growing terms in the trace formula |15| which do not permit 
to use standard methods. 
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The main purpose of this paper is to compute analytically the value of 
the 2-point correlation form factor, ^(t), in the limit r — > for two exam- 
ples of pseudo-integrable systems, (i) a plain billiard in the shape of right 
triangle with one angle equal ir/n and (ii) a rectangular billiards with a 
Bohm-Aharonov flux line inside. We argue that in the small-r limit the di- 
agonal approximation |16| is valid and the problem reduces to the calculation 
of distributions of periodic orbit lengths and areas occupied by periodic orbit 
families. Though for general pseudo-integrable systems very little is known 
on this subject, triangular billiards in the shape of right triangles with an- 
gle n/n belong to the so-called Veech polygons fI7|| , fl"8l and have a hidden 



group structure which make possible explicit calculation of necessary quanti- 
ties. After the calculations we found a finite value of the 2-point correlation 
form factor at the origin, < K 2 (0) < 1, which is different from both the 
Poisson distribution (for which K^Q) = 1) and the random matrix results 
(where K 2 (0) = 0). Analogous result has also been obtained for rectangular 
billiards with a Bohm-Aharonov flux line. Non-zero values of the 2-point 
correlation form factor at the origin confirm peculiar properties of spectral 
statistics for pseudo-integrable systems. We also discuss the comparison of 
theoretical predictions with the results of extensive numerical calculations. 

The plan of the paper is the following. In Section || the discussions of 
the trace formula and the diagonal approximation for the 2-point correlation 
form factor are presented. A brief introduction to the Veech structure of 
certain pseudo-integrable billiards is given in Section |||. For clarity we start 
in Section 3A_ with a simple example of square billiards where ideas and 
methods can easily be illustrated. Needed properties of the modular group 
and the Eisenstein series are shortly revised in Sections |3.2| and |3l| In 
Section [3.4| the Veech group for tt/u right triangle is derived and in Section 
3.4.2; the density of periodic orbits for this triangle is computed. In Section 
3.5 the calculation of 2-point form factor at the origin is performed and the 
comparison with the results of numerical calculations is discussed. Section 
|] is dwelt on the calculation of the 2-point form factor for a rectangular 
billiard with a flux line. As in the previous Sections the main point is the 
calculation of areas swept by periodic orbits around the flux line. The result 
depends on arithmetical properties of ratios of coordinates of the flux line 
to the corresponding rectangular sides. In Section |5| concluding remarks are 
presented. 
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2 The form factor in the diagonal approxi- 
mation 



2.1 The density of states 

The modern semiclassical approximation of multi-dimensional quantum sys- 
tems is based on various types of trace formulas which express quantum 
density of states (and other quantities as well) through quantities computed 
in pure classical mechanics |§, 0, [|Kj. The main step in deriving trace for- 
mulas is the semiclassical approximation for the (advanced) Green function 

G + (x,y) = ^ E _ En + ze , (3) 

n 

where E n is the set of energy levels and \l/ n the eigenf unctions as a sum over 
classical trajectories with energy E connecting initial point x and final point 

G+(x, y) = A tr exp(^S d - i-v). (4) 

tr 

S c i is the classical action computed along a trajectory, A tr is a pre-f actor 
depending on the system considered, and v is a phase (the Maslov index) 
which, roughly speaking, counts points where simple semiclassical approxi- 
mation breaks down. 

For 2-dimensional free motion (and for 2-dimensional polygonal billiards) 
the semi-classical approximation for G reads (see e.g. |J) 



where l p is the geometrical length of the orbit and k = v E is the wave vector 
(in the units h — 1 and m = 1/2). 

The knowledge of the Green function permits to find other quantum quan- 
tities as well. In particular the quantum density of states 

d(E) = Y,$(E - E n ) (6) 
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may be written by the means of the advanced Green function as 

d{E) = ~Im J dx G + (x,x). (7) 

The contribution from very short trajectories gives the mean level density, d, 
and the integration over the space selects periodic orbit contributions |J, || 
and determines an oscillating part of level density, d^ osc \E). For example, 
the density of states of an integrable rectangular billiard with sides a and b 
is 

d{E)=d + d {osc \E). (8) 

Here the smooth part is 

where A is the area of the rectangle (this formula is valid for all 2-dimensional 
billiards) and the oscillating part is 



oc 



p.p.o. n=l V V 



where 



l p = v / (2Ma) 2 + (2iV6) 2 . (11) 



In the rectangular billiard, the lengths of periodic orbits are 4 times degen- 
erate in the sum @ because (±M, ±N) give the same length. When the 
integral (|7]) is performed, orbits (M, N) and (— M, N) are absorbed in the 
same A p . The summation in (|T0|) is therefore performed over all primitive 
periodic orbits of length l p with M > repeated n times (an orbit (M, N) 
and its time-reverse companion (— M, —N) are counted as two different or- 
bits). In all integrable billiards periodic orbits are not isolated but belong to 
families. A p is the area of the pencil of periodic orbits of length l p . For the 
rectangular billiard A p = 2 A. 

Pseudo-integrable systems considered in the paper belong to the class of 
diffractive systems whose characteristic property is the existence of singular- 
ities which make the classical motion undetermined. Each time a classical 
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trajectory hits a singularity there is no unique way to continue it. Quantum 
mechanics smoothes out these singularities and associates with each (not too 
strong) singularity a diffraction coefficient, D(H,H'), (or scattering ampli- 
tude) which defines an amplitude of scattering on this singularity from the 
initial direction H to the final direction n' . 

Correspondingly, the semiclassical approximation of the Green function 
in the presence of a singularity at point x takes the form 

G(x, y) = G (x, y) + ^2 (^°> n)) D (ni n')G ((x , H'),y), (12) 



where Go(x,y) is the Green function without singularity and Go(x, (xo,n)) 
is a contribution to the Green function from a classical trajectory starting 
at point x and ending at the singularity xq with momentum in the direction 
n. Go((xo,n'),x') is a contribution to the Green function from a classical 
trajectory starting at point x with momentum in the direction n' and ending 
at point x! . 

This modification of the Green function changes the trace formula. For 
diffractive systems the density of states can now be written as the sum of 



three terms |g|, pTJ 

d(E) =d + d p . .(E) + d d . .(E), (13) 

where d is the mean level density, d p , Q , is the contribution of periodic orbits 
without singularity, and the third term, dd. .(E), is a contribution from all 
classical orbits starting and ending at the singularity (with, in general, dif- 
ferent momenta). These trajectories are called diffractive orbits and d± .(E) 
is a sum over all possible combinations of them 



d d .o.( E ) = y~] G(n 1 ,n' 1 )D(n l ,n 2 ) . . . G(n m _i, n' m )D(n' m , Hi) 

L — ' 7cm oE L — ' 



(14) 



where G(H,H') is the contribution to the Green function from a classical 
trajectory starting at the singular point with initial momenta in direction n 
and ending at it with final momentum in direction n' . 

For polygonal billiards the vertices with 7^ 1 play the role of singular 
point |19[. In the case of scattering on the angle a the diffraction coefficient 
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can be derived from Sommerfeld's exact solution fllS 



D(9 f ,9 i ) = - sin - 

7 7 



cos7r/7 — cos(6f + 9i)/ f y cos 71-/7 ~~ cos(#/ — 9i)/ r y 



(15) 



where 7 = a/n and Of (resp. 9j) is the final (resp. initial) scattering angle. 

For rectangular billiards with Aharonov-Bohm flux lines the flux lines 
themselves are singular points and the exact solution for an infinite plane 



with a flux line carrying a flux a |22j gives 



D(9 f A) = ^S-e^f-^. (16) 

COS -^j — 

The main difference between pseudo-integrable models discussed in this pa- 
per and usual diffractive models is the divergence of diffraction coefficients 
(|15D and (pl|) at certain directions (called optical boundaries because in the 
simplest case they separate illuminated regions from dark ones). Of course, 
exact solutions do not diverge even in vicinity of optical boundaries. The di- 
vergence comes from artificial separation of exact waves into geometrical and 
diffraction parts. Nevertheless, this formal divergence has profound effects 
on the structure of the trace formula. First, multiple diffraction along optical 
boundaries need a special treatment. Using a kind of uniform approximation 
]15| it was demonstrated that for polygonal billiards such multiple diffrac- 
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tion produces terms proportional up to a numerical factor to l/k where I is 
the total length of the diffractive orbit. When I is fixed and k — ► 00 (as in 
the usual approach to trace formulas) these terms are smaller than periodic 
orbit terms (pi]) but bigger than diffractive terms ([14]). But to compute spec- 
tral correlation functions one needs to consider a limit when k is fixed and 
I — > 00. In this limit multiple diffraction terms are bigger than both periodic 
orbit and diffraction terms. Another difficulty is related with the existence 
of terms corresponding to diffraction not exactly on optical boundaries but 
sufficiently close to them so their contributions are also large. Without exact 
summation of these quickly growing terms it is not possible to find spectral 
statistics of the systems considered. 

It the next Section we argue that, nevertheless, these terms give negligible 
contribution to the value of the 2-point correlation form factor at the origin 
and only diagonal contributions of periodic orbits will be important for this 
quantity. 
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2.2 The 2-point correlation form factor 

The 2-point correlation function is related with the level density by the formal 
expression 

Me) = (d(E+~)d(E-~)), (17) 

where the brackets denote an energy averaging around E on an energy win- 
dow much larger than the mean level spacing 1/d, and much smaller than 
energy E. 

The two-point correlation form factor is the Fourier transform of -R2(e) : 

/°° df i f f \ 

7 (d(E + -)d(E--))e 2 ^, (18) 

(the factors are chosen so that r and K 2 are dimensionless) . 

Trace formulas, roughly speaking, state that the density of states can be 
represented as a sum over classical orbits (both periodic and diffractive) 

d (osc) ^ = J- C p e iS ^ E)/h + c.c. (19) 
p 

Substituting this formal expansion into ([17]) and using the expansion 

S(E + e) « S(E) + T(E)e 



where T(E) is the period of classical motion one obtains ]16 



d(E + d(E - "-)) = (20) 
E C A < ex P l jt S vAE) - S P2 (E)) > e^+WW 

Here the terms corresponding to the sum of actions are omitted as it is 
assumed that they are washed out by the smoothing procedure. 

The corresponding expression for the 2-point correlation form factor is 
the following: 

K 2 ( T ) = £ ~c pi c; 2 < e «s n ms P2 (E))/H > s( t p1 + t P2 _ 2Mr ^ (21) 

P1,P2 
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The main difficulty in such an approach is the computation of the mean value 
of terms with action differences 



F(E) =< e i(s *'i (B) ~' s *'2(- B ))A > _ (22) 

The best developed approximation (called the diagonal approximation) con- 
sists in taking into account only terms with exactly the same actions |L6| 
i.e. 

pvzn _ / !> if S Pl (E) = S P2 (E) , , 

b W-\ 0, if S P1 (E) ? S P2 (E) ' ^ 

since terms with S P1 (E) ^ S P2 (E) will vanish by smoothing over E. In this 
approximation (assuming that for orbits with equal actions pre-factors are 
also equal (which is not always the case)) the 2-point correlation form factor 
takes the form 

Kt*\r) = J2 ^9 2 P \C P \ 2 5(T P - 2nhdr), (24) 
v 

where g p is the multiplicity of a given periodic orbit (i.e. the number of orbits 
with exactly the same action) and the summation is performed over orbits 
with different actions. In particular for integrable and pseudo-integrable 
systems from Eq. flIQ|) one gets 

Kf a9 \r) = ~J2 lj T^9 2 p S(l p - Ankdr), (25) 
p p 

where as before l p is the length of a periodic orbit and A p is the surface 
occupied by a periodic orbit family. 

It is instructive to perform the calculation for the simplest example of 
the rectangular billiard with sizes a and b. A periodic orbit in this billiard is 
defined by 2 integers m, n and its length is 

l p = ^(2ma) 2 + (2nb) 2 . (26) 

As pairs (m, n) and (m, —n) belong to the same family (or torus) the degen- 
eracy is g p = 2 (we remind the reader that in the rectangular billiard the 
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terms corresponding to m < are already taken into account in Ap), and it 
is sufficient to compute the density of periodic orbits with positive m, n 

p(l)= £ 8(1 -l p ). (27) 

m,n>0 

Changing the summation over integers (m, n) to the integration and using 
the substitution m = rcos0/(2a) and n = r sin 0/(26), one obtains by inte- 
grating over (j) from to 7r/2 

Pd) = u . (28) 

Since all families of periodic orbits in the rectangle cover the same area 
A p = 2 A and the length multiplicity is g p = 2, the 2-point correlation form 
factor for the rectangular billiard in the diagonal approximation is 

2A 2 f°° 1 

K 2 (t) = / J S ( l - Aixkdr ) P( l ) dl = !» ( 29 ) 

7T (I Jq I 

which is the expected value for the form factor of integrable systems [jTfifl. 
The diagonal approximation ( p3[) is known (with physical accuracy) to be 
valid for generic integrable systems |TB[ and can be modified P31 to compute 



mean values of more than 2 actions in the exponent of (p3|). 

For general systems the validity of the diagonal approximation is re- 
stricted only to small values of r |TfJ, and it is usually used to compute 



the first non-zero term of the expansion of the 2-point correlation form factor 
in powers of r. 

For diffractive systems with finite diffraction coefficient one can use the 
diagonal approximation for both periodic orbit terms and diffractive terms. 
But when the diffraction coefficient diverges in certain directions these cal- 
culations lead to difficulties. For example, multiple diffraction on optical 
boundaries corresponding to n repetitions of a primitive periodic orbit in 



pseudo-integrable billiard gives the following terms [15 



dmuit.di f f. (E) = ^2r c n cos(knl), (30) 



Ln 
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where c n are certain numerical coefficients. The attempt to use the diagonal 
approximation for these terms leads to the following result 

K (rnult.diff.) _ A ,2 r 3 ) (3^ 

if we take into account that the density of primitive periodic orbits in pseudo- 
integrable systems (at least for Veech systems (see the next Section)) differs 
only by a numerical factor from Eq. (f28|). But this expression contains pow- 
ers of momentum k and when k — > oo it cannot be correct. All terms corre- 
sponding to diffraction on or close to optical boundaries give similar quickly 
growing terms which cannot be treated separately. Without a ressumation 
of these terms the determination of spectral statistics of such models seems 
not possible. These arguments suggest the following scenario. The 2-point 
form factor is a sum of two terms 

K 2 (r) = f 1 (k a r) + f 2 (r), (32) 

where a is a certain positive quantity. The first function, fi(x), describes a 
result of resummation of quickly growing terms connected with divergence 
of the diffraction coefficient and when x — > oo fi(x) should go quickly to 
zero. The second function, fz(x), is a contribution of diffraction far from 
optical boundaries and can be computed similarly to ordinary diffraction 
in perturbation series of r. Of course, this is only a plausible conjecture 
and more detailed investigation should be done to give credit to it. 

Though the divergence of the diffraction coefficient prevents the calcu- 
lation of the 2-point correlation form factor in the full range, one can still 
use the trace formula (|14D to find its behavior at the origin, r = 0. The 
main point is that, even when the diffraction coefficient formally diverges 
the exact waves remain finite and using a uniform approximation (15 
can demonstrate that the ratio 



one 



D ^ n ' ] (33) 



kl 

is bounded for all angles, lengths and momenta. Each term in the diffractive 
trace formula (|l4]) is a product of certain number of these ratios and the 
total period of the corresponding composite orbit which appears due to the 
derivative over energy. Therefore it is of order of r multiplied by a constant 
and in the limit r — > all diffractive terms disappear. Only the periodic orbit 
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contribution (|T0| ) remains important at small values of r. From Eq. one 
concludes that for pseudo-integrable systems 

K2{0) = J^o g^J ^ 4^ 5(Zp ~ A7lkdT) - (34) 

The main problem now is how to compute the density of periodic orbits and 
the distribution of the areas of periodic orbit families. For generic pseudo- 
integrable systems very little is known and no reliable calculations can be 
done. E.g. for general plane polygonal billiards with angles commensurable 
with 7r it has only been proved |E7|, PBJ that the number of periodic orbits 
with length less than I, N(l p < I) obeys inequalities 

Cl l 2 < N{l p <l)< c 2 l 2 (35) 

for certain constants c\ and c 2 (depending on the polygon). But even the 
existence of an asymptotic law for J\f(l p < I) was not proved. 

Fortunately, there is a sub-class of pseudo-integrable billiards for which 
all necessary quantities can be computed due to the existence of a hidden 
group structure, and the triangular billiard in the shape of the right triangle 
with one angle equal to n/n belongs to this class. In the following Section 
we focus on those polygons. 

3 Veech structures for polygonal billiards 

We start the discussion of a hidden group structure of certain polygonal 
billiards with the simple example of the square billiard where the necessary 
ideas and methods can be illustrated clearly without technical difficulties. 

3.1 A simple case: the square billiard 

How can one evaluate the number of periodic orbits with length less than 
/ in a square billiard of size 2a with periodic boundary conditions? The 
exact expression for the length of the periodic orbits in the such billiard is, 
of course, 

l p = ^/(2ma) 2 + (2na) 2 (36) 
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with m G N and n 6 Z (see QITD). The number of periodic orbits with length 
less than /, M(l p < I), reads 

M{l p <l) =J2® ( l ~ 2aVm 2 + n^j (37) 



m,n 



and asymptotically when I — > oo 



/*oo /*oo . 

M{l p < I) = / dm dnQU- 2aVm 2 + 

JO J-oo ^ 



n 2 



,2 



< 38 > 

if one sets m = (rcos(p)/2a and n = (r simp) /2a. This is the number of all 
periodic orbits. More interesting questions and rich mathematical structure 
appear when one is interested in the calculation of the number of primitive 
periodic orbits M pp {l p < I) (that is, orbits with m and n coprime). 

The number of such orbits for a square billiard can easily be computed 
by using the inclusion-exclusion principle. The number of primitive periodic 
orbits with length less than I is the total number of periodic orbits with 
length less than I minus the number of orbits repeated p times with prime 
p, to which we add orbits repeated p±p2 times, which had been subtracted 
twice, etc. Finally one concludes that 

M pp {l p <l) = AT(/ p </)-^Ar(/ p <-) + ^AT(/ p < ' 



V — P\P2 

Pl,P2 F1F2 



E ^ r^r) • • • ( 39 ) 



P1P2P3 

Pl,P2,P3 



Using the I 2 dependence of Af in fl3"8|), we have 
^(1, < I) = M{l p <l){\ v ' ' w ' ^ ' 



where 



Z-e p2 (P1P2) 2 (P1P2P3) 2 

P ^ PUP2 ^ F1FZJ P1,P2,P3 ^^ l ^ A > 



CM = E ~, = If T^-r. < 41 > 
^— ' n s ±J - 1 — p s 

n=l p 
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is the Riemann zeta function. 
From (|38|) one gets 

V 2 

Kv(lv < = (42) 

Our aim is to generalize the previous calculation of Af pp (l p < I) to certain 
triangular billiards. This generalization naturally appears |l7j when one 
considers carefully the usual geometrical picture of the free motion inside the 
square billiard. It is well known that any trajectory of such a motion can be 
unfolded to a straight line when instead of the square billiard one considers 
the motion on the covering space which for square billiard is a plane with 
infinite square lattice of the side 2a. The vertices of this lattice (which are 
the images of the vertices of the initial square) have coordinates 

x = 2am, y = 2an (43) 

with integers m and n and can be considered as the result of the application 
of a 2 x 2 matrix with integer coefficient to a horizontal vector (2a, 0) 

/ m k \ ( 2a \ / 2am \ , . 

U I J { ) = {2an)- ^ 

Thus, the periodic orbit lengths ( |36D are the distances between these vertices 
and the initial point (0,0). 

The problem of finding the number M pp of primitive periodic orbits with 
length less than / is therefore equivalent to the problem of finding out how 
many 2x2 matrices with integer coefficients and determinant equal to 1 (since 
m and n are coprime one can impose ml — nk = 1) exist with m 2 + n 2 < x 2 
for a given x (or, which is equivalent, with n 2 + 1 2 < x 2 ). The group of 2 x 2 
matrices with integer coefficients and determinant equal to 1 form a group 
SX(2,Z) and in the next two Sections we shall discuss its main properties. 
Though this material is well known we find it useful to remind it informally. 

3.2 The modular group 

The subgroup of SX(2,R) containing all 2 x 2 matrices with integer coeffi- 
cients and determinant equal to 1 is called the modular group SL(2, Z). The 
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standard representation of this group (see e.g. P6"| ) is the Poincare half-plane 
Ti. with measure 

ds 2 = ^(dx 2 + dy 2 ) ; (45) 

y2 

A matrix g G SL(2, Z) is represented by the isometry 

g : H —> H 

mz + k (46) 



z 



nz + I 



The modular group is generated by the translation T : z i— > z + 1 and the 
inversion 5 : z i— > — l/z, which correspond respectively to the matrices 



1 a 
1 



(with a = 1 for the modular group) and 

1 
-1 



(47) 



(48) 



Since the modular group is a discrete group, we can define its fundamental 
domain T> (shown in fig. [I]), that is the domain of the Poincare half-plane TC 
that covers H, under the action of the representation ( fiq) of the group. 

In order to compute the number of matrices g = ( m ^ J G SL(2, Z) 



verifying n + I < x , we have to evaluate 



AT(x)= Yl 1 ( 49 ) 

g e roc\G 

n 2 + l 2 < x 2 

where G = SL{2, Z) and is the subgroup of G generated by the trans- 
lations (r^ = {T n ,n G Z}): since the left multiplication by matrices of the 
form T p 

1 p \ f m k\ ( m + pn k + pi \ , . 

1 Jin I = n Z (50) 
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-1 



-1/2 



1/2 1 X 



Figure 1: The fundamental domain of the modular group 

does not change n and I it is necessary to consider the quotient T^G (i.e. 
2 matrices which differ by T v are considered only once), so that the sum 



is convergent ||29|| . If we assume that in the limit x — > oo the sums can be 
written as integrals over n and I with uniform measure (see later) in the form 
(B/n)dndl, we get 

J\f(x) = [ -dndl = %x 2 . (51) 

3.3 Eisenstein series 

In order to compute the coefficient B, let us introduce the Eisenstein series 

E(z,s)= Yl ( Im 9(z)) s (52) 
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for s > 1. From expression (pE6| ) we get, since ml — nk = 1 



1m g(z) = V —- (53) 

\nz + 

where y = Im z. Since Im g'g{z) = 1m g(z) for </ G Too, the sum over 
roo\G is well defined. Let us first compute the asymptotic behavior of E(z, s) 
when s — > 1. For a given R G R, we can rewrite the sum ( |52"D as a finite sum 
over elements of G for which n 2 + I 2 < R 2 and a sum over elements for which 
n 2 + I 2 > R 2 which diverges as s — > 1. The divergent part, n 2 + I 2 > R 2 , is 

£ diV (M = - f (54) 
TT J n 2+i 2 >R 2 \nz + l\ 2s 

By s f°° f n r l - 2s drd6 . , 

(55) 



7r Jq [(x sin + cos 0) 2 + (y sin 0) 2 ] s 
Since 

oo E>2(l-s) i 

r i-2 Srfr = " „ (56) 



2(8-1) «-i 2(s-l 



and the finite part of the Eisenstein series can be neglected as compared with 
the divergent part, we have 

E (z s) ~ yB ^ (57) 

»~i 2tt{s-1)J (xsin0 + cos0) 2 + (ysin0) 2 ' 1 j 

The computation of the integral can be performed the following way : setting 
A = 1 — x 2 — y 2 , B = 2x and C = 1 + x 2 + y 2 , we get 

2 



(x 2 + y 2 ) sin + cos 2 + 2x sin cos Jo A cos 20 + B sin 20 + C 

(58) 

and since C 2 — A 2 — B 2 = Ay 2 , the integral is equal to 

/ ; = - (59) 

Jo VA 2 + B 2 cos^ + C y 



Finally 



E(z,s) ~ — (60) 

v ' ' 2(s-l) K ' 
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and this limit does not depend on z. 

Now let us integrate the series in (]52f) with the invariant measure d/j,(z) = 
dxdy/y 2 over a part, V Y , of the fundamental domain V corresponding to a 
restriction y <Y. If o? is the width of the fundamental domain, 

/ E(z,s)dii(z) = ^2 (Im g(z)) s dfi(z) 
Jvy ger ao \G J ' Dy 

= £ / (Xm z) s dpt(z) (61) 
dx\ Y %y s (62) 

Jo y 

d- — -. (63) 



In transformation from Eq. fl6"I|) to Eq. (|62|) we take into account that the 
image of fundamental region T> (and T>y) under the action of G is a certain 
region on the upper-half plane which under the action of can be moved 
into a vertical strip of width d (which is the fundamental region for the group 
Too). These images can not intersect and when Y — > oo will cover the whole 
strip with y <Y. 

The asymptotic behavior of this integral is thus the following 

lim / E(z, s) dfi(z) = . (64) 

s-* 1 Jv Y 8 — 1 

By comparing this expression with Eq. (B^) one concludes that the value of 
the constant B is 

B = 2^— (65) 
Vol V y 1 

and from Eq. ([51]), the final answer for the density of primitive (n and / are 

coprime) periodic orbits in a square billiard is 

K P (n 2 + l 2 < x 2 ) = ^^x 2 (66) 

For the modular group d — 1, Vol T> — 7r/3, and for a square billiard with 
side 2a we have x = I /(2a) according to ( |36l) , so 

3/2 

KvHv < l) = ^ (67) 
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which agree with Eq. ( pE2| ) obtained by a different method. 



3.4 Veech group for ixjn right triangles 
3.4.1 The symmetry group 

The previous calculations were possible because we have found a group -the 
modular group- that relates periodic orbits in a square to a simple vector 
(see Eq. (fH|)). In order to generalize this construction for more complicated 
polygons it is important to point out that the modular group is the symmetry 
group of the unfolding of the square billiard (that is, the lattice whose unit 
cell is a 1 x 1 square). Indeed this square lattice is, evidently, invariant under 
the following two transformations: the rotation by 7r/2 around the center of 
the square which we denote by S and the translation of one coordinate (say x) 
by 1 which we denote by T. In Cartesian coordinates these transformations 
are represented by the following matrices 

*-Ui) T -(iD- (68) 

The vertices of the lattice are also unchanged under the action of the group 
generated by S and T. But it is well known that this group is exactly the 
modular group SX(2,Z). Therefore the modular group plays a double role 
for a square billiard. First, it is the group of invariance of unfolding of the 
square and, second, it generates periodic orbits starting from a fix vector as 
in Eq. ©. 

It has been proved by Veech [17] that for certain polygons (called the 
Veech polygons) there exists a group with similar properties. In particular, a 
n/n right triangle (i.e. a triangle with angles ir/2, ir/n, n{n — 2)/2n) belongs 



to the Veech polygons pTfl , [18 



Let us consider this case in details. The geometrical construction of the 
unfolding of the classical trajectories in such a billiard is slightly different for 
n even and odd. By reflections with respect to the sides corresponding to 
the it /n angle the n/n triangle can be unfolded to the regular n-gon. For 
n even the opposite sides of this n-gon should be identified by translations 
(see Fig. 0a ). For n odd one has to consider 2 regular n-gons reflected with 
respect to one side and to identify parallel sides by translations as in Fig. |2|b. 
The resulting surface is the surface of genus (n — l)/2 for n odd (see 0) to 
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a) 



b) 



Figure 2: The unfolding of n/n triangle. Left - n is even. Right - n is odd 

which all trajectories belong. From this construction it is clear that if a 
group of invariance exists it should include the rotation by 2ir/n around the 
center of these n-gons. In Cartesian coordinates this rotation is defined by 
the following matrix 



To find other transformations which leave this surface invariant it is necessary 
to consider a few families of periodic orbits. 

For n even, we define in the n— gon two important elementary families 
of periodic orbits: the first one is the family of horizontal primitive periodic 
orbits, the second one is the family of primitive periodic orbits making an 
angle n/n with the horizontal (see fig. |3|). For n odd we only define the first 
family (see fig. £|). 

For n even (n = 4p or n = 4p + 2), the first family has orbits with lengths 




(69) 




n 



n 



(70) 
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and widths 

^■ = 2cos^^sin^ (71) 
n n 

with 1 < j < p. The second family has orbits with lengths 

(2j - 2)tt 7i 

I' = 4 cos ^ '— cos - (72) 

J n n 

and widths 

, (2j - 2)tT . 7T 

w,- = 2 cos sin — (73) 

J n n 

with 2 < j < p if n = 4p or 2 < j < p + 1 if n = Ap + 2. The orbit with j = 1 
is special: it has a length and a width equal to 

l[ = 2cos7r/n, = 2shi7r/n (74) 

For n odd (n — 2p + 1), the lengths and widths are the following 

t a ■ 2 3 n 71 (K\ 
Lj = 4 sin cos — (75) 

n n 

and 

W j = 2 sin — sin-. (76) 

n n 

with 1 < j < p. It is of interest to compute the ratio of the length of each 
periodic orbit family to its width. From the above formulas it follows that 
for all families, except the one with j = 1 for even n, this ratio is the same 

- = 2cot-. (77) 

w n 



For the exceptional family (74) this ratio is 2 times smaller. 

The unfolding of any family of periodic orbits gives an infinite strip of 
points of period L p and width W p . If there is a group of invariance of the 
unfolded surface, it should include a transformation which leaves invariant 
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Figure 3: The elementary orbits in the decagon 



periodic orbit strips. Assume that the strip is oriented horizontally. In this 
case one sees that the shift of the form 

{li) < 78 > 

leaves points of the strip invariant provided that 

L p = naW p , (79) 

where n is an integer. 
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Figure 4: The elementary orbits in the heptagon 



Because for periodic orbits considered the ratio (|77f) is constant the in- 
variance group should include the following transformation 

( 1 2cot^ \ /on , 
r„= Q (80) 



Veech proved JET], |18| that the invariance group for 7c/n triangle is a dis- 



crete subgroup T n of SL(2,WL) generated by the two elements (|69D and 
Similarly to the relation ([0]) for the square, periodic orbits in this triangle 
are generated by the action of T n over the elementary families of periodic 
orbits considered above. 

We shall call the members of these families "basis orbits" (and label them 
by the index i). We define the corresponding basis vectors Vi by Vj = (Lj, 0) 
and = (ljCOs(ir/n),l'jSm(ir/n)), so that {v^ = {vj,v'A for n even and 
{vi} = {vj} for n odd. 

An element g of the symmetry group T n has the following matrix repre- 
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sentation 



The result of the action of this element to one of the basis vectors Vi = 
(vUfVii) gives the coordinates of a new primitive periodic orbit (more pre- 
cisely, a periodic orbit situated on the boundary of periodic orbit pencil) 

^ = (82) 



c d J \ v i2 

and the length of this primitive periodic orbit is the length of this vector. 

The first family of basis vectors for even n and for all basis vectors for odd 
n can be chosen in the form = (Li,0) and the lengths of periodic orbits 
generated by applying the group T n is 



L g = Va 2 + c 2 Li, (83) 



where a, c are matrix elements of g fl8"2"|). The second basis periodic orbits, 
v'i, are obtained from horizontal vectors by rotation by ir/n. But the matrix 
corresponding to the inverse of this rotation 



cos - sin - 



sin - cos - 

n n 



(84) 



does not belong to our group T n . Nevertheless, this matrix plays the role of 
a Hecke operator, namely, even if it does not belong to T n the conjugation of 
any matrix from this group does belong to T n : if g G T n , then r~ x gr G T n . 
To prove it let us note that 

r-V£r = al (85) 

where a n is the generator (|69|) because all rotations commute, and it is easy 
to check that 

r _1 r n r = -o n T~ x , r~ X T~ l = T n <J~ x . (86) 

The right-hand sides of these relations belong to T n and as all matrices from 
T ra can be written as a product of generators we get r~ l gr G T n for g G T n . 
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Using this conjugation one can rotate the second family of periodic orbits 
for even n by —ix/n and the lengths of the orbits generated by the vectors v\ 
will be related to matrix elements of T n by the same relation as in flS"3| ) 

L' g = Vtf + fl'i. (87) 

Therefore, to find the distribution of periodic orbit lengths it is necessary 
to compute the distribution of a 2 + b 2 for matrices from T n , which has been 
done for the modular group in the previous section. Eq. (^) can be derived 
the same way for T n . According to the previous section one can compute the 
density of periodic orbits and other quantities as well by investigation of the 
fundamental domain of T n . 

The distribution of areas of periodic orbit families is also easy to obtain: 
as all matrices from T n have unit determinant, the area covered by the pencil 
corresponding to g V{ (g G T n ) is the same as the area covered by the pencil 
corresponding to basis vectors Vi, i.e. it is equal to LjWj. In other words, 
there is a one to one correspondence between pencils of primitive periodic 
orbits and vectors g Vi for g e T n . The discrete group T n is related to 
periodic orbits in the n/n right triangle in the same way as the modular 
group is related to periodic orbits in the square. 



3.4.2 The density of periodic orbits 

The fundamental domains of the symmetry groups T n for n even and odd 
are described in Figs. [5] and |] respectively. For n even, it is the union of two 
triangles with angles 2ir/n, 0, on the Poincare half-plane : the area of the 
domain is Vol T> = 1-n(n — 2)/n, and its width is 2cot(7r/n). 

For n odd the two triangles have angles ir/n, tt/2 and 0, therefore Vol T> = 
ir(n — 2)/n; the width of the fundamental domain is 2 cot n/n. These shapes 
of fundamental domains can be obtained by taking into account that the 
group T n considered as a group acting on the Poincare upper-half plane as 



in Section |3]2] includes (i) the translation by 2 cot w/n and (ii) the rotation 
around point i by 2it/n for n odd and by 4ir/n for n even. This difference 
between even and odd n is related to the fact that the rotation by angle ir cor- 
responds to the transformation g \— > —g, but these 2 matrices are represented 
by the same function on TC (see ([46])). For odd n (n = 2q + 1), the group 
generated by the generator (|69|) contains rotations by angles 27rj/(2g+ 1) for 
j — 0, 1, . . . 2q. The value j — q + 1 corresponds to the rotation by is + n/n. 
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Figure 5: The fundamental domain of T n for n even 



As the rotation by tt is the identity transformation, the rotation by ir/n be- 
longs to T n and it is a primitive generator of the subgroup {cr^,p G Z} of 

rv 

For even n (n = 2q), the rotation by 2irq/n is the identity, therefore the 
rotation by ir/n does not belong to the group and the primitive generator of 
the subgroup {cr^,p G Z} is the rotation by 2n/n. 

Due to (|66|) , we now get the number of matrices 



a b 
c d 



G T T 



(88) 



verifying a 2 + c 2 < x 2 



M{a 2 + c 2 < 



x 



n 



71 



cot — X 
7T(n — 2) n 



2n 



TT 



cot — X 
k n(n — 2) n 



(n even 



(n odd) 



(89) 
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Figure 6: The fundamental domain of T n for n odd 

These formulas give the total number of matrices verifying a 2 + c 2 < x 2 . 
But due to the existence of rotation matrices in the group T n each primitive 
periodic orbit length appears a few times in the above calculations. This 
multiplicity corresponds to different unfoldings of a given periodic orbit. 

For n odd the 2n matrices of the form {±f3^g, < k < 2n — 1}, where 
(5 n is the matrix of rotation by 7r/n, give rise to one primitive periodic orbit. 
For n even there exist n matrices of the form {±a^g, < k < n — 1} (where 
a n is the matrix of rotation by 27r/n (^)) which describe one periodic orbit. 

For g given by (|8"ED, the length of gvi is L^a 2 + c 2 . So the number of 
primitive periodic orbits of type gv^ less than / is 

^ iL ' < i) = ^Yr'HrS ■ m 

The number of all primitive periodic orbits is the sum over all such con- 
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tributions: 

M pp {L p <l)=C 

where 

A= .... 

4 n 



M PP (L P <l)= C-, (91) 



1 • 27r 

A=-sm— (92) 



is the of our triangle and 

e4 < 93 ) 



I n 7T 

COS 



27r(n-2) n 
For n odd (n = 2p + 1) Eq. (|75|) gives 

^ If ~ 16 cos 2 n/n ^ sm 2 (2nk/n) ' ^ 

For n even (n = Ap + 2e, epsilon = 0, 1) from Eqs. (0), (0), and (|74|) it 
follows 



1 _ 1 l» 1 ^ I 

^ L 2 16cos 2 7rM \ ^-f cos 2 (2? - lWra ^ 



16cos 2 7r/n \~[ cos 2 (2j — l)n/n ^—^ cos 2 {2j — 2)n / n 

'(n-2)/2 



+ 4 



E 27 / + 4 • ( 95 ) 



16cos 2 7r/n I ^-^ cos 2 kix/n 

The last sums can be calculated using the evident formulas (for another 
method of calculation see [|17]]) 

1^1 

sin x ' (x — 07r)^ 



and 



1 1 xV I 1 



7T 2 ^^2o-l-x) 2 + (2o-l + a;) 2 ^ ^ 



cos 2 7ra;/2 n 2 ^ K (2q - 1 - x) 2 {2q - 1 + 
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Taking into account that sin 2 (£;7r/n) = sin 2 ((n — k)ir/n) for odd n and per- 
forming the following transformations 



n— 1 co ^ ^ ^ 2 | 

fc=l q=— oo v ' t— — oo q=—oo 

one obtains that for odd n 

El rr — 1 

sm 2 (2nk/n) = 6 ' ^ 



Similarly for even n 



(n-2)/2 2 

1 n z 



Therefore the value of constant C is 

C = r v I K ~ol f r ° dd " • (101) 

1927r(n - 2) \ (n 2 + 20) for even n v ; 

This result corresponds to primitive periodic orbits with geometrically dif- 
ferent lengths: time-reversal orbits are not included in the summation. They 
give an additional factor of 2 in Eq. ( |101| ). In [T7| the orbits corresponding 



to different unfoldings of the same periodic orbit have been included in the 
asymptotic formula, which leads asyptotically to the additional factor n for 
even n and 2n for odd n in Eq. ( |101|) . Furthermore, if one needs all peri- 



odic orbits including repetitions Eq. (|101|) should be multiplied by n 2 /6 as 
in Section pTT| . 

In Fig. [7| we present numerical results of the cumulative density of primi- 
tive periodic orbits in it/ 8 right triangular billiard (with area A = 4n) when 
all orbits (time-reversal and for different unfoldings) are included. The solid 
line is the best quadratic fit to these data 

M PP (L P <l) = .0294/ 2 - .6055/ + .3617. (102) 

One sees that this fit can hardy be distinguished from numerical results. 
The theoretical prediction for the coefficient in front of I 2 is C/A according 
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to Eq. (|9"I|) ; C is given by Eq. ( |101| ), and has to be multiplied by n since the 
numerical computation has taken into account the repetitions of each orbit, 
and by 2 as time-reversal orbits are taken into account as well: 



C_ 
A 



24vr 2 



.0295, 



(103) 



which is in excellent agreement with numerical calculations. 
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Figure 7: The cumulative density of periodic orbits for the 7r/8 right trian- 
gular billiard. 



3.5 Explicit calculation of the 2-point form factor for 
the ir/n right triangle 
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3.5.1 First case: no degeneracy of the lengths 

We assume in this section that there is no degeneracy between the lengths of 
the periodic orbits (except the ones connected by the time-reversal transfor- 
mation) or more carefully, that there is no pair of primitive periodic orbits 
whose lengths are commensurable. Since the lengths of the gvi are propor- 
tional to the lengths of the Vi, the necessary requirement for the validity of 
this condition is the absence of commensurability relations between the Lj. 
In this case the 2-point correlation form factor in the diagonal approxi- 



mation is done by Eq. (|25|) . The sum over different periodic orbits can be 
split into a sum over all types of periodic orbits, then the sum over periodic 
orbits of each type can be replaced by an integral; since the density pi of 
periodic orbits of type i only takes into account once the periodic orbit and 
its time-reverse, the degeneracy is g p = 2 and 

POO 

Y.9l = A H dl p^) • ( 104 ) 

p i J o 

In (|25|), A p is the area occupied by a pencil of periodic orbits of length L p ; 
but this area is the same for all trajectories belonging to the same family i. 
So we just have to evaluate the area Ai = LiWi occupied by an elementary 
orbit of type i. The lengths fl70|), fl72|) and flf|) are 

kit it 

Lfc = 4 cos — cos — , 1 < k < p — 1 (105) 
n n 

and L = 2 cos7r/n if n = 2p, and 

2klT 7T 

L fc = 4sin cos—, l<k<p (106) 

n n 

if n — 2p + 1. The widths are only half the widths Wi given by (|71~1), ([73]) 
and (|76| ) since each fundamental pencil is symmetric with respect to the line 
joining two images of the 7r/2 corner of the triangle (see fig. |3] and ||). So 

the Wk are 

/l7T 7T 

W k = cos — sin-, < A; < p — 1 (107) 

n n 



32 



if n = 2p, and 



2 /l7T 7T 

W k = sin sin-, 1 < k < p (108) 

n n 



if n = 2p + 1 . Then for small r 

1 2 /-ex) 



/12 /-oo i 

K 2 {r) =4V-i -= 5 (Z - Ankdr) Pi (l)dl, (109) 

o7T Jq la 



and replacing the density of orbits of type % by its mean value 

_ 7T 2 dAf iiPP (L p < I) 
Pl ~ 6 dZ 



(110) 



where Ni tPP (L p < I) is given by Eq. fl9"U| ) we obtain, when performing the 
integral, 

where the average density of states is J = .4./47T with 

1 2i , . 

,4= -sin — 112 

An y ' 

is the of the triangle. 

The sum (|109|) over the widths fllOTQ and (|108|) gives 

f^sin 2 -, n even 
4 n ' 

(113) 



So we finally get 



with 



! ^sin 2 -, n odd 

4 n ' 



*.(r) = -, + e{n { (H4) 



e(n) = when n is odd, (H5) 
e(n) = 2 when n is even. (H6) 
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3.5.2 Second case: degeneracy of the lengths 

We have assumed in the previous Section that the lengths of all primitive 
periodic orbits were non-commensurable. In the case of the n/n right tri- 
angle, there may exist a commensurability relation between the Lf. given by 
( |105|) or ( |106| ) if there is one between the cos{kn/n) (0 < k < p — 1) for n 



even, or between the sin(2fc7r/n) (1 < k < p) for n odd. It is shown in 
that if n is an odd prime, there is no such relation between the sin(2A;7r/ri). 
Ref. (|]31|) deals with the case (k,n) = 1 and gives the same conclusion. It 
seems that in the general case, only one relation of that kind exists between 
the cos(/c7r/n), which is 

71 7T 

2cos(--) =cos(0) (117) 
3 n 

and that no relation exists between terms with sinus. Therefore the only 
degeneracy occurs in the case where n is even and 3|n, that is n 6 6Z. In 
that case, from fllCf) we get 

d{E + |) d{E - ^)) = K dia9 + 

y AA^ 1 i H i p -i pl )+i%(i p +i pl ) + (ng) 

p+^p + v 

where K dtag is the usual diagonal approximation (^1|). p + means that we 
only count for one orbit in the sum the orbit and its time-reverse, therefore 
there is a coefficient 4. 

If there is a relation rriiLi = m 2 L 2 (with mi and m 2 coprime) between 
two lengths of primitive periodic orbits, we have a contribution Rf g to the 
2-point correlation function ([[7]) which comes from orbits of lengths qm\L\ 
and gm2-^2, g G Z* : 



2 A-ttZ 9-ttJc ^ 



2 



= 2AlA2 V 1 e ^ (gmiLl) (119) 
47r 2 .27rA; ^ orniLx ' v ; 

The sum over all repetition numbers q of a function of gLx (where is a 
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primitive periodic orbit) can be replaced by an integral 



R 



dog 



AlA 2 

4iT 3 k 



dl- 



mil 



(120) 



where the density pi(/) is the density of periodic orbits of type 1 (that is, with 
length proportional to L{) with length less than I, given by ( |110| ). Performing 
the Fourier transform (|T8|) and the integral over I gives 



K 2 (0) 



cot 7r/n 
67r{n - 2)d 



w x w 2 

Wim 2 



(121) 



In our case the degeneracy is given by (|117|) and 



We finally obtain 



W Wn = I sin 2 -. 
— 2 n 
3 



(122) 



K 2 (0) 



n + e(n) 
3(n-2) 



(123) 



with 



when n is odd 
e(n) = \ 2 when n is even and 3 f n 
6 when n is even and 3 I n. 



(124) 



This formula is the main result of our calculations for the triangular billiards. 
It clearly demonstrated the peculiarities of spectral statistics for pseudo- 
integrable systems. The non-zero value of the form factor (< 1) at the origin 
does not correspond to any random matrix ensemble but it is typical for 
intermediate statistics [25R, p2|. 



3.6 Comparison with numerical calculations 

To compare the prediction (|123| ) with numerical results we have computed 
20000 levels for triangular billiards in the shape of a right triangle with one 
angle n/n for all n — 5, 7 . . . ,30 (the case of n = 6 is integrable). For each 
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Figure 8: The 2-point form factor for 7r/8 right triangular billiard. 



triangle we take levels from 15000 till 20000 and compute numerically the 
corresponding 2-point correlation form factor. A typical result is presented 
in Fig. ^. From data like this, it is quite difficult to find the value of the form 
factor at the origin because r — > corresponds, according to Eq. (|18"1), to an 
infinitely large energy difference in the 2-point correlation function: therefore 
numerically we always have K 2 (0) = 0. We found it convenient first to fit 
the numerical data to the following simple expression for the form factor, 

K 2(T) = *-^f;\ (125) 



and then from it compute fG(0). 



K 2 {0) = 1 - -. (126) 
a 
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The form ( |125| ) has been chosen because (i) one wants a simple expression, 
(ii) when r — > oo the form factor should go to 1, (iii) to describe the level 
repulsion it is necessary that 

OO 1 

(1 - K 2 (r)) dr = (127) 
and (iv) the expression ( |125| ) when a = 4 equals the form factor of the 



so-called semi-Poisson model p5" |, p2| which serves as a reference point for 
intermediate statistics. 

We stress that the above expression has no solid theoretical explanations 
and it is used because it relatively well describes our numerical results. The 
only fitting parameter is -fG(O) related with a by Eq. (|126|) . We tried two 
fitting procedures. First we fit Eq. ( |125| ) for the data with all r or, second, 



to decrease the influence of very small r, where numerical accuracy is not 
very good, we did not consider the data with < r < 0.25. In Fig. || these 
two fits are presented. The first one gives 1^(0) ~ 0.44 and the second one 
K 2 (0) « 0.565. The expected value ( p3| ) for n = 8 is 5/9 « 0.56. 

In Fig. ^ the results of such fitting procedures are given for all triangles. 
Lower two curves correspond to these fits and the upper curve is the predic- 
tions ( |123|) . (Of course, only points are important. Curves are presented for 
clarity.) The numerical results quite well follow theoretical formula (|123|) but 
there is a small shift which decreases when the region of small r is ignored. 
This difference between the curves seems to be a consequence of the fact 
that the result ( |123| ) corresponds to asymptotic limit k — ► 00 but numer- 
ical calculations have been performed at large but finite energy. To check 
this point we present in Fig. [XOj the results of the calculation of the mean 
number variance, E^(L), for the 7r/30 right triangle, in 10 energy intervals 
[8000/c, 8000(A: + 1)], < k < 9 (the energy increases from bottom to top). 
It is well known (see e.g. [|IJ]) that the behavior of the mean number variance 
at large distances is related with the value of the form factor at the origin 
by simple formula 

S (2) (L) -> K 2 (0)L when L -> 00. (128) 



From Fig. |10| it is clearly seen that even for 80000 levels the curve does not 
stabilize. To find its limiting behavior we extrapolate point by point (with 
L fixed) this ten curves with a fit A(L) + B(L)/^/k. (It means that for 
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Figure 9: K 2 (0) for n/n right triangles, n = 5 to 30. Circles are theoretical 
results ( |123|) . Squares are the fit ( |125| ) when the region of small r, < r < 



.25, is omitted. Diamonds are the same fit but with all r. 



each L we fit 10 points to find the best A(L) and B(L).) The limit curve 
(i.e. A(L)) is the most upper curve in Fig. [H| It perfectly reproduces the 
expected features of E^ 2 ^(L): it is a straight line with slope -^2(0) = 0.38 
corresponding to the expected value (|123| ) for n = 30. In the same way, 
Prosen and Casati []33| have computed TP\L) for triangle billiards with 
angle ir/5 for much larger values of the energy, and it seems that such fit 
works well for their calculations and the result for ^2(0) agrees with ( |123| ). 
These (and other) calculations clearly demonstrate that the value of the 2- 
point correlation form factor at the origin converges slowly to the theoretical 
result ( |123| ) with increasing energy. This behavior may be a consequence of 
the conjectured existence of two different terms ( [32]) in the form factor and, 
in the final extent, a manifestation of the strong diffraction in vicinity of 
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Figure 10: 'E^(L) for energy windows with higher and higher energy 



optical boundaries. 



4 Rectangular billiard with a flux line 

4.1 Preliminary calculations 

This section is devoted to the study of a rectangular billiard with the Aha- 



ronov-Bohm flux line \22\ at a point r~o = (xo, yo) inside the rectangle. In the 
polar coordinates, r, ip, around this point the vector potential of the flux line 
has only ip component 



A, 



a 



(129) 
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and the 2-dimensional Schrodinger equation for the motion in this potential 
is (when h = c = 1 and m = 1/2) 



— + -— + — (— -iaV -f / ' 
dr 2 r dr r 2 \dtp 



* n (r,(p) = 0. (130) 



Similarly to triangular billiards discussed in previous Sections this model 
belongs to the class of diffractive systems. The diffraction coefficient for 
the scattering on the flux line (|1~6|) diverges in the forward direction but as 
for pseudo-integrable systems the contribution of diffractive orbits can be 
neglected when computing the value of the 2-point correlation form factor at 
the origin. 

It is well known that the Aharonov-Bohm potential ( |129| ) does not change 
classical trajectories but gives an additional phase, A0, when a trajectory 
turns n times around the flux line 

A0 = 2mra. (131) 

Therefore the contribution of periodic orbit to the trace formula (|T0|) will con- 
tain an additional phase depending on the winding number of the trajectory 
around the flux line. 

Periodic orbits in the rectangle of sides a, b are determined by two integers 
M and N in the usual way and they are characterized by their length 

l p = ^{2Ma) 2 + (2Nb) 2 , (132) 

the area occupied by the periodic orbit family, and the winding number 
around the flux line. Each pencil of primitive periodic orbits occupies an 
area 2ab = 2A, so its width is 2A/l p . The images of the flux line in the 
unfolding of the rectangular billiard are located at the points 

((( 1 + 2k)a,(( 2 + 2k')b). (133) 

Here Q takes the values q or 2 — (i = 1, 2), where 

ei = — e 2 = — (134) 
a b 

are the ratios of coordinates of the flux line to the corresponding sides, and 



k, k' G Z (see Fig. |TT|). Let us define [x] as the largest integer less than or 
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Flux line 



Periodic orbit (M=l, N=2) 



The same orbit unfolded 



Pencil of periodic orbits 



Figure 11: An unfolded trajectory in the rectangle 



equal to x, so that 



[x] < x < [x] + 1, 



(135) 



and {x} = x — [x] G [0, 1[. 

Each unfolded pencil of a primitive periodic orbit contains two and only 
two images of the flux line, since it covers twice the area of the rectangle. 
The periodic orbits from this pencil parallel to the vector (M, N) and going 
through the images of the flux line (which we shall call the saddle connections) 
split the pencil of primitive periodic orbits parallel to (M, N) into three 
pencils of same length (see Fig. [TTJ) . Only the central strip is affected by the 
presence of the flux line and any trajectory from this strip gets a phase 2ira 
(according to (I3~l) and since the orbit is primitive). So the winding number 
of a periodic trajectory is nothing but the repetition number of a periodic 
orbit belonging to the central pencil. 

To compute the widths of the central strip as a function of M and N, let 
us note that the algebraic distance from an image + 2k) a, (C2 + 2k')b) of 
the flux line to the saddle-connection linking the points (0, 0) to (2Ma, 2Nb) 
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IS 

2A 

d=— (C2M - CiiV + 2k' M - 2kN) ; (136) 

Lp 

The two images of the flux line which are inside the pencil are the two nearest 
points among the all images (((i + 2k)a, (C2 + 2fc')6) with integer k and k! that 
have positive distance to the saddle- connect ion. They correspond to points 
such that the distance fll36f) is positive and less than 2A/l p , or less than 1 
in units of 2A/l p . 
Let us set 

Q± = [e 2 M ± eiN], e± = {e 2 M ± e x N}. (137) 

The four possible values of (Ci, C2) ((61,62), (61, 2 — e 2 ), (2 — 61,62), and 
(2 — ei, 2 — 62)) give four possible families for the distance ( |136| ) : 

di = 2ki±{Q + + t + ) 

<U = 2ki±(Q- + eJ) (138) 

where ki is a certain integer which depends on k, k', M and N. Among 
these four families, exactly two points correspond to a distance positive and 
less than 1 : for instance if Q + and Q_ are even, only 2k\ + Q + + e + and 
2/c 2 +Q_+e_ can be positive and less than one. So we must have 2ki+Q + = 
and 2&2 + Q- = and the two images of the flux line that are in the pencil 
of primitive periodic orbits are at a distance e + and e_ from the saddle- 
connection (0, 0) — (2Ma, 2Nb) ; if both Q + and Q_ are odd, the two distances 
are 1 — e + and 1 — e_ . Since the width of the central strip is the difference 
between the two distances, it is in both cases |ej — 62 1. Dealing in the same 
way with the case where Q + and Q- have opposite parity, we get that the 
width of the central strip in units of 2A/l p is 

e_ — e + | , if Q + and Q_ have the same parity 

|1 — e_ — e+| , if Q + and Q- have opposite parity 
Both cases can be summed up in the following formula : 

T} = tp(x,y), (140) 
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with 

x = e 2 M + €\N 

y = e 2 M-e 1 N ' { 1 

where 

y) = \f(x) - f(y)\ , f(x) = (-l)W ({x} - §) (142) 

f(x) is an even function of period 2 ; if we restrict the study of ip to [—1, 1] x 
[—1,1] we have 

, >. f \x + y\ if xv < /+ in \ 

<p(x,y) = < , ., y ~ _ , (143) 

^ v ' y; I x-j/ ifij/>0 ' v ; 



and the Fourier expansion of ip is 



(cos 7rnx — cos irny) 
V (x,y) = 2j2~ — 2 — • 144 



2 



n=l 



Using ( |141|) and ( |144j) , we obtain that the width of the central strip associated 
to the orbit (2M, 2N) is the following 

8 sin 2 (7me2M) sin 2 (7meiiV) .-\ 

• ( 145 ) 

n=l 

4.2 Form factor for the billiard with flux line 

The value of the two-point correlation form factor in the diagonal approxima- 



tion given by Eq. (25), when diffractive contributions have been neglected, 
still holds for billiards with flux line. But A p now includes the phase fac- 
tor depending on the repetition number of the trajectory. The density of 
periodic orbits fllCf ) becomes 

1 00 A 1 

= 5 E Et^-^^ (146) 

pp+,pp~ n=l V P 

Here we distinguish between two types of orbits. The orbits associated with 
a primitive orbit pp + have a phase exp(2i7rno;) for the orbit repeated n times, 
and the total coefficient in the trace formula associated with these orbits is 
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where A Pl , A P2 and A P3 are the areas covered respectively by the three strips 
in which the pencil of periodic orbits splits. The orbits associated with 
a primitive orbit pp~ have the complex conjugate phase exp(— 2iirna) and 
their contribution to the trace formula is proportional to 

•Ap- n A Pl ~\~ Ap 2 c + Ap 3 . (148) 

When the terms with same length l p = nl pp are gathered together, we get 

d P AE) = £ ^-l=e^-^-^ + c.c. (149) 

where 

A np = A Pl + A P2 cos(27rna) + A P3 (150) 



and the sum ^2 p + goes over orbits M, N > 0. Equation ([H]) now becomes 

pp + n=l yl x ' 

If rj p is the width of the central strip expressed in units of 2A/l p , we have 

A np = 2A(l — i] p + T] p cos2iina) = 2A(1 — 2r) p sm 2 7ma), (152) 

and (using the fact that d = A/An) the 2-point correlation form factor at 
small t is a sum of 3 terms 

8A^ 1 ^ 1 f Akr 

n=l PP PP v 

32.4 sin 2 irna y—^ i] pp f Akr 



E sm una rj pp I 

n 2 ^ Irm V 

n=l pp P p v 

32^4 *^ sin 4 Ttna v-^ ^™ r ( i Akr 



+ -^r 2^— ^--rr ■ (153) 



/ v 9 / v ( \ VV 

7r z — ' n z L — ' l vv \ n 

n=l pp ^ l 

The summation over primitive periodic orbits can be done by replacing the 
sum by an integral, taking into account the density of primitive periodic 
orbits. If e\ and €2 are rational numbers 

e l = (154) 
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where pi and qi are co-prime integers, the width of the central strip ( |145|) 
only depends on the remainder t\ of M modulo q\ and r 2 of N modulo q2 

8 °° 1 

^(ri,r 2 ) = — > — sin 2 (7r — nri) sin 2 (tt — nr 2 ). (155) 

There are q\q<i periodic orbit families 

M = qi k + r 1} N = q 2 k' + r 2} (156) 

with k, k! e N. 

To compute sums in Eq. (|153|) one needs to know the mean density of 
primitive periodic orbits for each family, Ppp, ri ,r 2 (l) • 

Let c be the greatest common divisor of q±, q 2 . c = (q%, q 2 ). If (r\, r 2 , c) 7^ 
1, then M and N are not coprime and there is no primitive periodic orbit. 
In the opposite case it is demonstrated in Appendix A that 

Ppp,n,r 2 {l) = Pp P (l)a(r 1 ,r 2 ), (157) 

where p pp (l) is the mean density for all primitive periodic orbits in the rect- 
angle (cf. Eq. (0)) with M, N > 

V 

MO = 4^. (158) 

and 

<*{ri,r 2 ) = = 1 J] (1--) J] (1--), 

(159) 

and lcf(gi, g 2 ) is the least common factor of qi, q 2 . 

The knowledge of the mean density of periodic orbit families permits the 
computation of mean values of different quantities depending on families. If 
f(ri, r 2 ) is such a quantity its mean value is defined as follows 

</>= J2 /(n,r 2 )«(n,r 2 ). (160) 

Ti mod qi 

(r 1 ,r 2 ,c) = l 
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In particular 



E 



I I rrn 



Akr 



pp 



n 



^2 V /3 ( r i,r 2 ) 



Ti mod qi 
(r 1 ,r 2 ,c) = 1 



(1)5 I 



Ah 



n 



A-kA 



(161) 



< rf > . 



The sums over n that appear in ( |153| ) can be computed using the standard 
formula 



cos 2imx 



^ (277//) 

n=l V 1 



X 



1 

X+ «' 

6 



for < x < 1. 



(162) 



It gives 



oo . o 9 

sin 7ma 7r _ . _ N 
_ = T a(l-a) 

n=l 



E 



(163) 



and 



oo . 4 

sin ?rna 



E 

n=l 



7T" 



-a: 



(164) 



where a is the fractional part {a} of the flux through the rectangle when 
< {a} < 1/2 and a = 1 - {a} when 1/2 < {a} < 1. 

Using ( |161|) , (|163|) and ( |164j ) one concludes that the 2-point correlation 
form factor for r — ► ( |153| ) is the following 



K 2 (0) = 1 - 12a(l - a) < r] > +6a < rj 2 > 



(165) 



To use this formula it is necessary to know the values of < i] > and < rf >. 
In the case where both ei and e 2 are irrational non-commensurable quantities 
the fractional parts {ne 2 M} and {neiN} cover the whole interval [—1, 1] and 
rj and rf can be computed by integrating expression ( |143|) of rj over the square 
[— 1, 1] x [—1, 1]. Simple calculations show that in this case 



1 3' 



(166) 
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< V 2 >= \- (167) 

Therefore when the coordinates of the flux line are non-commensurable with 
the corresponding sides 

K 2 (0) = 1 - 3a + 4a 2 . (168) 
In Appendix B, it is shown that for all rational q 

< V >= \, (169) 

like in the irrational case. The average < rj 2 > is more difficult to compute 
analytically: we have found an analytical expression only when q\ divides q 2 
(or similarly q 2 divides q\). 

Though the general formula for < rj 2 > is cumbersome, the computation 
of the < rf > for rational e\ and e 2 can easily be done numerically using 
Eqs. (|141|) and (|143|) . For small denominators the values of < rj 2 > are given 



in Table [IJ. 

To check the obtained formulas we have computed numerically 1500 first 
energy levels for the rectangular billiard with sides a = 4 and b = n and 
the flux line with coordinates (from low left corner) x = 5a/9 and t/q = 
116/20. The typical picture of K 2 (t) is shown in Fig. |12[ As for the triangular 
billiards discussed in previous Sections we extrapolated K 2 (r) to small r 
using the simple expression ( 125|) . The results for different values of the flux 



are presented in Fig. [T3|. We also check the following more suitable fit (which 
obeys the condition ( 127 ) when c = (1 — b) 2 ) 



K 2 (r) = \ b Y T ^ hen T n (1 ^ )/C ■ (170) 
K ' \ 1, when r > (1 - b)/c K ' 

It gives practically the same results. The existing numerical precision does 
not permit to distinguish these 2 fits. 

In the case where xq = 5a/9 and yo = 116/20, simple calculations show 
that < rf >= 4867/29160; Eq. Qgjg ) gives the expected value of K 2 (0) 



K 2 (0) = l-^^-a + 4a 2 . (171) 
v ; 4860 v ; 
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7 


8 


2 


1/3 


2/9 


1/4 


2/9 


2/9 


2/9 


11/48 


3 


2/9 


2/9 


13/72 


17/90 


5/27 


47/252 


107/576 


4 


1/4 


13/72 


1/6 


61/360 


1/6 


85/504 


1/6 


5 


2/9 


17/90 


61/360 


14/75 


89/540 


37/210 


167/960 


6 


2/9 


5/27 


1/6 


89/540 


4/27 


35/216 


47/288 


7 


2/9 


47/252 


85/504 


37/210 


35/216 


26/147 


85/504 


8 


11/48 


107/576 


1/6 


167/960 


47/288 


85/504 


1/6 



Table 1: Value of < rf > for a rational flux line 
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Figure 12: The 2-point form factor for the rectangular billiard with flux line 
with a = 0.4 and its smoothed value (white line). 

which corresponds to the solid curve in fig. [13] (Note that the coefficient 
of a equals approximately 2.99 and is practically indistinguishable from the 
coefficient 3 for irrational e»). 

Similarly as for triangular billiards there is a small difference between 
the theoretical and numerical curves. For triangular billiards where more 
levels are available this difference slowly decreases with energy. We expect 
the same behavior also for rectangular billiards with a flux line. 

5 Conclusion 

In this paper we have obtained explicit expressions of the 2-point correlation 
form factor K 2 (t) in the limit r — > for a few typical examples of pseudo- 
integrable billiards : triangular billiards in the shape of right triangles with 
one angle equals ir/n, and rectangular billiards with a flux line. The obtained 
values of K 2 (0) differ from standard examples of spectral statistics (the Ran- 
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Flux 

Figure 13: ^(0) for different values of the flux a (points) and the asymptotic 
theoretical prediction (solid line). 

dom matrix theory and the Poisson statistics), which confirm analytically 
the peculiarities of spectral statistics of pseudo-integrable systems. The cal- 
culations have been performed by analysing analytically the properties of 
classical periodic orbits of the systems considered. 

In order to elucidate further the special properties of spectral statistics 
of polygonal billiards, it would be of interest to compute ^(0) for generic 
triangular billiards without the Veech structure. Moreover, we have taken 
into account only the diagonal terms and, consequently, were able to obtain 
only ^(0). The computation of the next terms in the expansion of if 2 Or) in 
powers of r should include the exact ressumation of singular contributions, 
coming from the diffraction close to the optical boundaries. The solutions 
of these problems require the development of new methods beyond the ones 
used in this paper. 
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Appendix A 



Periodic orbits in a rectangle with sides a, b are determined by 2 integers M 
and iV which count the difference of coordinates of initial, (xi,yi), and final, 
(xf,Vf), points 

Xf = Xi + 2aM, y f = y l + 2bN. (172) 
The length of the periodic orbit is the geometrical length of this vector 



L p = v/(2aM) 2 + (2bN) 2 . (173) 

The mean cumulative density and the corresponding quantity for primitive 
periodic orbits (when M, N are co-prime integers) can be computed as for 
the square billiard (see Eqs. (|38|) and (f42"D ). When I — > oo and if only positive 
M are considered 

N(L r < - ^, (174) 

and 

V 2 

N ^ L - < ss- (175) 

The purpose of this Appendix is the computation of the mean cumulative 
density of primitive periodic orbits for periodic orbit families when 

M = 7*1 mod qi, N = r 2 mod q^. (176) 

The asymptotics of N pp (L p < I) when / — > oo is related with the behavior at 
small x of the B-function associated with these periodic orbits 

e(x) = ^e- x ^. (177) 

pp 

If 

C 

Q{x) -> — , when Z -> 0, (178) 
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then 

C 

N PP (L P < 1) -> ^^y z27 ' when z 00 • ( 179 ) 
We are interested in the following 0-function 

oo 

6( X )= ^ e - 3 ,((2aM)»+(2W0») j (lg0) 
M,N=-oo 

where the summation is performed over all integers M, N with the following 
constraints 

(M, N) = 1, M = n mod g 1? N = r 2 mod g 2 - (181) 

Note that both positive and negative values of M, N are considered. When 
only positive M are taken into account the formulas below have asymptoti- 
cally factor 1/2. 

To impose the restriction M = r mod q it is convenient to introduce the 
6-f unction 

( n' iH 7,° • (182) 

,y 10, otherwise 

Its explicit form may be the following 

1 9-1 

5 1 = - y e 2mkt/q . (183) 
q / — 1 

As in Section RTTj the condition (M, iV) = 1 can be taken into account by 
the inclusion-exclusion principle 

oo oo 

f(M,N)= Yl ^2f(Mt,Nt)ti(t), (184) 

(M,N)=1 M,N=-oo t=l 

where fi(t) is the Mobius function equal (—1)™ if t is a product of n distinct 
primes, if t contains a squared factor, and = 1. 
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Combining all the necessary restrictions one finds the final expression for 
the 6-function (ff5D|) 



oo oo 

0(x ) = _l j2 E j2^ e ~ AxtHM2a2+N2b2) 

M,N=-ooki mod t=l 
xe 2mk 1 (Mt-r 1 )/q 1 +2Trik 2 (Nt-r 2 )/q 2 (185) 



Using the Poisson summation formula 

oo _ 
^ e -xM 2 +2niyM = I* ^ ^(M+yf/x^ ^ggj 

one obtains that 



. x 

M=— oo M=— oo 



^ /fW p -27rifciri/gi-27rifc 2 r2/g2 



Anh'Tn-.n^ E/ ' 

fci mod t=l 

x x 

X 

M=-oo N=-oo 



e -^ 2 (Af+fcit/ Q i) 2 /(4xt 2 a 2 ) e -i T \N+k 2 t/ q2 f /(4xt 2 b 2 ) ^ lg7 ^ 



When x — > the dominant contribution comes from terms with zero expo- 
nent, i.e. from terms with 

M + ^ = 0, andiV + ^ = 0, (188) 

or kit = mod q{. The asymptotics of the 0-function is therefore the fol- 
lowing 

Q ( x ) = Tir— E E^^W^e- 2mfciri ^~ 2 ^ 2 ^. (189) 
Aabx qiq 2 , ^— f-r' ^ 

fei mod qi 5=1 

Using the representation (|183|) for these (^-functions and performing the sum- 
mation over ki one gets that when x — > 

©(^Ht^Ws), ( 19 °) 
Aabx 
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where 



1 

F(r 1 ,r 2 ) = V V<5 /lt _ riigi ^ 2t _ r2)(?2 -— . (191) 

Q1Q2 , A , 7 1 

li mod qi t=l 

From Eq. (|1 79|) one concludes that the asymptotics of the mean cumulative 
density of primitive periodic family ( |181| ) (with M, N > 0, i.e. with a factor 
1/4) is 

N pp (L p <l) = ^- b F( ri ,r 2 ). (192) 

To perform the summation over t in Eq. ( |191| ) it is necessary to know the 
number of solutions of two equations 

lit = ri mod qi, l 2 t = r 2 mod q 2 . (193) 

It is well known (and can be easily checked) that the number of solutions of 
the equation ax = b mod g depends on the greatest common divisor of a 
and q, (a, q) = d. If d = 1 there is 1 solution, x = 6a _1 mod g. If d > 1 
and c? f b there is no solutions. If d\b there is one solution, x = (b/d)(a/d)^ 1 
mod (g/cf) and consequently, there are d solutions modulo q: Xj = x +(q/d)j , 
j = 0, . . . d — 1. Therefore 



F(r 1 ,r 2 ) = ^^^-{qut){q 2 ,t)5 {qu t), ri S( q2 ,t) 



M2% 



1 



t 2 



'■■2 



rfi|(<?i,fi) (gi,t)=di 

d 2l(92- r 2) (92.*)= d 2 

Terms corresponding to (dj, da) > 1 give a contribution to the sum, since in 
that case q%, q 2 , d\ and d 2 have a common factor, which contradicts condition 



(M,N) = 1. The sum (|194|) can therefore be restricted to (di,d 2 ) = 1, and 
the sum over t is now a sum over i' where t = d\d 2 t' '. Let us denote by P the 
product of the prime factors of q\ that do not divide c and by P' the product 
of the prime factors of q 2 that do not divide c. Now 

(c,P) = (c,P') = (P,P') = l. (195) 
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If p is a prime dividing di and c = 52), then p divides di = (qi,t). Since 
it divides c it also divides 52, so p\(q2, t) and (di, d 2 ) 7^ 1, which is impossible. 
So the prime divisors of d\ have to be taken among the divisors of P, and in 
the same way the prime divisors of d 2 have to be taken among the divisors 
of P' . Similarly, we can check that if p is a prime divisor of c which divides 
t', p divides d\ = (qi,t), so d\ would contain a prime factor of c, which is 
impossible; and if p is a prime factor of P or P', p divides both d x and d 2 . 
So the sum over t' must be restricted to t' which do not contain any prime 
divisor of q± or q 2 . As fi(ab) = ji{a)jj,{b) for co-prime a and b, one gets 

F(r 1 ,r 2 ) = — ^ d x d 2 -—— 2 (196) 

qiq2 d M d ^ 

where the sum is taken over all d±,d 2 , t' verifying di\(P, r{), d 2 \(P', r 2 ), t' \ q±, 
t' \ q 2 . Using the identity 



n(1 1, 

p\k S\k 



we get 



Hn,r 2 ) = — n n ^--^ n ^--^ 

H H Pt<?l,pt92 F p\(P,ri) F P\(P',T2) F 

= J\(l-^)a(r u r 2 ) = - 2 a(r u r 2 ), (198) 

J. J. pZ n z 

all p 



where 



a(ri ' r2) = ^n 1 n-i/r^ II (!-^) n t 1 -^)' 

9l92ll p | lcf(gii(?2) U VP J p | ( , 1>ri)iptc P p|(, 2 ,, 2 )^c P 

(199) 

and \cf(qi,q 2 ) is the least common factor of q±, q 2 ()cf(qi,q 2 ) = Q'l^/c). 

It is also instructive to check directly that at{r\,r 2 ) are normalized cor- 
rectly 

a(n,r 2 ) = l, (200) 

n mod <ji 

(n ,t-2,c)=i 
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where as above c = (qi,q 2 )- We use once more the inclusion-exclusion prin- 
ciple 



E /(n,r 2 ) = $>w E ^> r2 )- ( 201 ) 



Ti mod qi t\c ri mod g; 

(''l,'-2,c) = l t|rj 



If we set 



D = qi q 2 J] C 1 "^)' ( 202 ) 



P 

p|lcf (<?i,<?2) 



we get, with ([L9|) and (gOlj) , 

E «(™> = ^EM«) E E E ^ E " ( ' 



D ^— ' ^ ^ ^ 5i ^ 5 2 

Ti mod gi t\c ri(t|ri) r 2 (i|r 2 ) <5i|(P,ri) <52 1 (J 3 ',^^) 

( r l.''2' c ) = 1 

= ^E/^)E E E E ^r^r 

t|c <5i|P <5 2 |P' ri(i<5i|n) r 2 (« 2 |r 2 ) 

= ^EirE^pE^|- ( 203 ) 

t\c <5i|P 1 5 2 |P' 2 

Here we have used the fact that if t\c and 5\P, t and 5 have no common 
factor. In the above sums it is always understood that the summation over 
Ti goes only from = to — 1. But the last sum in Eq. ( |203|) exactly 
equals D because cPP' = lcf(gi,g2) and Eq. ( P00| ) holds. 



Appendix B 

In the same way one can compute the mean value of r] defined in Eq. (115-51 ) 

8 °° 1 

<V>= E 7 l(n,r 2 )a(r 1 ,r 2 ) = —J2 — J2^ ( 204 ) 



n 2 D ^— f n 2 

n mod qi n=X t\c 

(f1,r a ,c)=l 



V* M^OM^) . 2 ( rA . 2 / r 2 
y > > > — ^ — sin I 7rn — | sin I 7m — 



<5l #2 V Ql J V ^2 

n(t\ ri ) r 2 (t|r 2 )<5i|(P ) ri)<5 2 |(P' ! r 2 ) \ ^/ \ ^. 
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Since 

^si n 2 7m- = -^(l-(W, ? ), (205) 



q 2t5 

tS\r 



one obtains 



t\c 5i\P 1 8 2 \P' 2 

oo . 

/] — i 1 - <W,«i) (1 - S ntS2m ) . (206) 

n=l 

The sum over n includes 4 terms. The first is the sum over all n 

00 i 2 

E^=V < 207 > 

n=l 

The second sum has the restriction that n = {qi/t5\)m and 

?1=1 

The third sum is the same but with the substitution 1 — > 2. The fourth 
sum incorporates two restrictions, nt#i = mod q\ and ntd^ = mod g2- 
Remembering the definition of P and P' (see (|195| )) one concludes that in 
this last case the restriction is n = (cPP'/(£5i5 2 )) m and 



n=l ^ ' 



Performing the summation over S{ and t in Eq. ( |206| ) one notes that all three 
last sums will have factor 

<5i|P 5 2 |P 



But for any K > 2 we have 



$>(*) = 0. (211) 

S\K 
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Since q\ and q2 are greater than 1, it is impossible that simultaneously c = 
P = 1, or c — P' — 1, the terms (|210|) equal zero. Therefore only the term 
( |207| ) survives and ( |206} ) gives 

z . 1 \- M*) \- M^i) \- M^a) , 01 _* 

< »7 >= LirLipLx ( 212 ) 

t|c <5i|P 1 5 2 \P' 2 

These sums are exactly equal to D and finally we get 

<V>=\- (213) 
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